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In this study, we use the cell dynamics method to test the validity of the Kormogorov-Johnson- 
Mehl-Avrami (KJMA) theory of phase transformation. This cell dynamics method is similar to the 
well-known phase-field model, but it is a more simple and efficient numerical method for studying 
various scenarios of phase transformation in a unified manner. We find that the cell dynamics method 
reproduces the time evolution of the volume fraction of the transformed phase predicted by the KJMA 
theory. Specifically, the cell dynamics simulation reproduces a double-logarithmic linear KJMA plot 
and confirms the integral Avrami exponents n predicted from the KJMA theory. Our study clearly 
demonstrates that the cell dynamics approach is not only useful for studying the pattern formation 
but also for simulating the most basic properties of phase transformation. 
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1. Introduction 

Phase transformation occurs by the nucleation and subsequent growth of a nucleus in a 
system where the first-order phase transformation takes place. It has attracted much attention 
for more than a half century 1-5 from a fundamental point of view as well as from technological 
interests. These include the mechanical properties of metallic materials, 6 the recrystallization 
of deformed metals, 7 and the manufacturing of basic thin-film transistor devices, such as solar 
cells and active matrix-addressed flat-panel displays. 8 

The nucleation and growth processes are often described in terms of the old standard the- 
ory called the KJMA theory developed by Kolmogorov, 1 Johnson and Mehl, 2 and Avrami. 3-5 
According to this theory, the time evolution of the volume fraction of a new transformed 
phase follows the linear KJMA plot with the integral Avrami exponent n that is given by the 
slope. However, it is recognized that this theory often fails to explain experimental results; 9, 10 
neither the KJMA plot becomes linear nor the Avrami exponent becomes an integer. There 
is also some debate about the validity of the assumption used in the theory. 11 To resolve the 
discrepancy, a realistic yet efficient simulation method that could take various factors into 
account is indispensable. 
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The direct atomic-scale computer simulation of the kinetics of phase transformation using 
the molecular dynamics or Monte Carlo method is still a difficult task. Even the most funda- 
mental phenomenon, such as nucleation, is not easy to study using these methods. Instead, the 
problem of phase transformation has been tackled using a coarse-grained Ginzburg-Landau- 
type model called the Cahn-Hilliard, 12 Ginzburg-Landau 13 or phase-field 14 ' 15 model, which 
requires the solution of highly nonlinear partial differential equations. Since this model re- 
quires the time integration of the partial differential equations, it is not easy to simulate the 
long-time behavior of the dynamics of phase transformation 16 except for special traveling wave 
solutions. 17,18 

A model based on a cellular automaton instead of a partial differential equation can im- 
prove the efficiency of numerical integration. 6, 19 It has been used to test the KJMA theory 
in the recrystallization of metals. Although the cellular automaton is sufficiently flexible to 
implement the various local reactions of recrystallization, it does not have a direct connec- 
tion to the equilibrium phase diagram. Therefore, the connection between the phase diagram 
and the phase transformation is not so clear compared to the Ginzburg-Landau-type model 
mentioned above. 

In this study, instead, we use the cell dynamics method 20,21 to study the validity of the 
KJMA theory. This method is attractive because it has the merit of cellular automata and 
is computationally efficient, and yet it keeps the connection to the phase diagram through 
the Landau-type free energy. The format of this paper is as follows: in §2, we review the 
cell dynamics method and present the necessary modification for studying the nucleation and 
growth. In §3, we follow closely the work by Jou and Lusk 14 and test the validity of the KJMA 
theory using the cell dynamics method. We conclude in section 4. 

2. Cell Dynamics Method for Nucleation and Growth 

To study the phase transformation, it is customary to study the partial differential equa- 
tion called the phase-field model 14, 15 which is equivalent to the time-dependent Ginzburg- 
Landau (TDGL) 13 or Cahn-Hilliard model: 12 

at [> 

where 5 denotes the functional differentiation, ip is the nonconserved order parameter, and T 
is the free energy functional. This free energy is usually written as the square-gradient form 

?W\ = \j [£(W) 2 + h{$)] dr. (2) 

The local part h(ip) of the free energy T determines the bulk phase diagram and the value 
of the order parameter in equilibrium phases. The double-well form was frequently used for 
h(tp) to express the two-phase coexistence and study the phase transformation between these 
two phases. 
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This TDGL equation (1) for the nonconserved order parameter ip was loosely transformed 
into a space-time discrete cell dynamics equation by Puri and Oono 21 following a similar 
transformation of the kinetic equation for the conserved order parameter called the Cahn- 
Hilliard-Cook equation. 20 In their cell dynamics method, the partial differential equation (1) 
is replaced by a finite difference equation in space and time in the form 

iP{t + l,n) = F[^(t,n)], (3) 

where the time t is discrete and an integer, and the space is also discrete and is expressed by 
the integral site index n. The mapping F is given by 

Fty(t, n)} = -f^{t, n)) + [< n) » -^{t, n)] , (4) 

where f{ip) = dh(ip) / dip , and the definition of <C • • • 3> for a two-dimensional square grid is 
given by 20-22 

j=nn i=nnn 

where "nn" denotes nearest neighbors and "nnn" next-nearest neighbors. An improved form 
of this mapping for a three-dimensional case was also obtained. 22 

Oono and Puri 20 ' 21 further approximated the derivative of the local free energy f(ip) called 
the "map function" in the tanh form 

/(V) = ^ ~ V-^4tanh^, (6) 



23 



with A = 1.3 that corresponds to the free energy 

h(ij>) = -A In (cosh V) + ^ 2 (7) 
and is the lowest order (0(ip 2 )) approximation to the double well form of the free energy 

K^) = -\^ + \^ (8) 

when A = 1.5. They 20 ' 21 used this simplification since this cell dynamics system is invented 
not to simulate the mathematical TDGL partial differential equation but to simulate and 
describe the behavior of nature directly. Later, Chakrabarti and Brown 23 discussed that this 
simplification is justified since the detailed form of the double-well potential h(ifj) is irrelevant 
to the long-time dynamics and the scaling exponent. 

Subsequently, however, several authors used the map function f{tp) directly obtained from 
the free energy h(tp) as it is 24 ' 25 and found that the cell dynamics equation (3) is still tractable 
numerically. Ren and Hamley 25 argued that one can easily include the effect of the asymmetry 
of the free energy and the asymmetric characteristic of two phases using the original form of 
the free energy function f(i/j). 

Despite the popularity of this cell dynamics method in the soft-condensed matter com- 
munity as a simulator of pattern formation due to various factors, it has not yet been used 



3/16 



Jpn. J. Appl. Phys. Regular Paper 

to study the most fundamental problem of phase transformation by nucleation and growth. 
In the next section, we use a parameterized free energy function, study the kinetics of phase 
transformation, and test the validity of the KJMA theory. 

3. Numerical Results 

3.1 Two-dimensional growth of single domain 

To simulate the growing stable phase after the nucleation, we have to prepare the system 
in a state where one phase is metastable and has the higher free energy than the other stable 
phase. The free energy difference between the stable and metastable phases is determined 
from the supersaturation in liquid condensation and from the undercooling in crystal nucle- 
ation. Microscopically, this free energy difference is necessary for the nucleus to overcome the 
additional curvature effect caused by the interfacial tension and to continue to grow. 26 

To study the growth of the stable phase using the cell dynamics method, we consider the 
time-dependent Ginzburg-Landau (TDGL) equation (1) with the square gradient free energy 
functional (2). The local part of the free energy h(ip) we used is 14 



4 IT v r/ 2 V 3 2 , 
This free energy is shown in Fig. 1, where one phase at tp = is metastable while another 
phase at tp = 1 is stable. The free energy difference Ah between the stable phase at tp = 1 
and the metastable phase at tp = is solely determined from the parameter e: 

Ah = h(ip = 0) - h{ip = 1) = ^. (10) 

Therefore, e represents the supersaturation or the undercooling. 

The metastable phase at ip = becomes unstable when r\ = 3e, which defines the spinodal. 
The height AE of the free energy barrier at tp = (rj — ?>e)/2n can be tuned by the parameters 
n and e: 

^-*(*- 2 ir)-«*-°>- 2L: ^- (11> 

which vanishes when n = 3e at the spinodal. 

The steady-state analytical solution of the TDGL with a constant interfacial velocity 
has been obtained in one dimension by Chan 17 when the free energy h{tp) is written in the 
quadratic form such as in eq. (9). Using Chan's formula, the interfacial velocity v of our TDGL 
model (1), (2) with the free energy (9) is given by 




(12) 

Chan 17 further suggested that if the interfacial width is small, the interfacial velocity of a two- 
dimensional circular or three-dimensional spherical growing nucleus is asymptotically given 
by eq. (12) of the one-dimensional model. The larger the free energy difference e and the lower 
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Fig. 1. Model double-well free energy (9) that can realize two-phase coexistence when e = 0. The 
parameter e determines the free energy difference Ah, and the parameters r\ and e determine the 
free energy barrier AE. 



the free energy barrier r\ are, the higher the front velocity v becomes. 

The critical radius R c of a two-dimensional circular nucleus is also given analytically: 14 ' 17 

R<^ = ^. (13) 

v 3e 

In the metastable phase, the circular nucleus of the stable phase with a radius (R) smaller than 
R c shrinks, while that with a radius larger than R c grows and its front velocity approaches 
eq. (12). Again, the larger the free energy difference e and the lower the free energy barrier r] 
are, the smaller the critical radius R c is. 

We implemented the above free energy (9) into the cell dynamics code written by Math- 
ematica TM 27 for the animation of spinodal decomposition developed by Gaylord and Nishi- 
date, 28 and simulated the growth of a single circular nucleus of a stable phase. 

Initially, we prepared a small circular nucleus of a stable phase within a metastable phase 
and simulated the growth of that nucleus. The system size is 100x100 cells, D in eq. (2) is 
D = 0.5, and the periodic boundary condition is imposed. The initial order parameter ip is 
randomly chosen from 0.9 < ip < 1.1 for the circular nucleus of the stable phase and from 
—0.1 < ip < 0.1 for the metastable environment. The diameter of the initial nucleus is fixed at 
d = 11. Therefore, the initial nucleus occupies a part of 11x11 cells. The effective area of the 
stable phase is computed by counting the number of cells with the order parameter ijj > 1/2. 

Figure 2 shows the effective radius of the circular nucleus of the stable phase calculated 
from the effective area of the nucleus as a function of time step t. The nearly linear growth 
of the nucleus of the stable phase is clearly visible, which indicates a constant front velocity 
for the stable-metastable interface predicted from the analytical solution for the TDGL. 17 
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Fig. 2. Evolution of effective radius of area of stable phase plotted as a function of time when rj = 1 
and e = 0.1 (open circles), n = 1.0 and e = 0.3 (open triangles), and r] = 0.3 and e = 0.1 (open 
diamonds). The effective radius shows a linear time dependence, indicating a constant interfacial 
velocity. 



The velocities v estimated from Fig. 2 and predicted from eq. (12) are compared in Table 
I together with the critical radius R c calculated from eq. (13). It can be seen that eq. (12) 
correctly predicts the general trend of the front velocity v when the two parameters r\ and e 
are altered. Since the cell dynamics method does not solve the TDGL directly, the discrepancy 
between cell dynamics simulation and theoretical prediction (12) by a factor of roughly 2 is 
not very significant. 



Table I. Comparison of front velocities v estimated by cell dynamics simulation and prediction using 
eq. (12) for various potential parameters r\ and e. The corresponding critical radii R c are also 
tabulated. 



T] 


1 


1 


0.7 


0.5 


0.4 


0.3 


e 


0.1 


0.3 


0.1 


0.1 


0.1 


0.1 


v (simulation) 


0.063 


0.22 


0.082 


0.10 


0.12 


0.14 


v [eq. (12)] 


0.15 


0.45 


0.18 


0.21 


0.24 


0.27 


Rc [eq. (13)] 


3.3 


1.1 


2.8 


2.4 


2.1 


1.8 



From the above comparison of the values estimated by cell dynamics simulation and 
theoretical prediction using the steady-state solution of the TDGL for a two-phase system, 
we consider that this cell dynamics method is effective for studying the phase transformation 
by the growth of the multiple nucleus of the stable phase. 
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3.2 KJMA kinetics by cell dynamics simulation 
3.2.1 Site-saturation nucleation 

In site-saturation nucleation, a fixed number of nuclei are prepared initially, and subse- 
quent growth is monitored. The KJMA theory gives an analytical expression for the volume 
fraction / of the stable phase as a function of time t. In two dimensions, the formula leads 
to 14 

/ = 1 - exp (-Trn v 2 (/ + t ) 2 ) , (14) 

where no is the number density (number per unit area) of the randomly distributed initial 
nuclei, v is the growth rate of the radius of each nucleus discussed in the previous section, to is 
the origin of time which can hopefully take the incubation time of nucleation into account. 15 
Prom eq. (14), we have 

log (- ln(l - /)) = 2 log (t + t ) + constant. (15) 

Therefore, the KJMA theory predicts that a double logarithms log (— ln(l — /)) versus 
log (t + to) is a straight line that is known as the KJMA plot with the integral tangent n = 2, 
which is called the "Avrami exponent" . 
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Fig. 3. Evolution of volume fraction / calculated by cell dynamics simulation for site-saturation 
nucleation as a function of time (open circles) when n = 1.0 and e = 0.3. The broken line denotes 
the theoretical prediction (14) with to = 0, while the solid line denotes that with to = 10. A better 
agreement between the simulation and theoretical results is attained when the incubation time 
to = 10 is taken into account. 

We have simulated the site-saturation nucleation using the cell dynamics method. Now 
a finite number of nuclei of the stable phase is distributed over the area we considered. The 



7/16 



Jpn. J. Appl. Phys. Regular Paper 

initial nuclei are circular and have the diameter d = 8, which is larger than the critical 
radius R c in Table I. Then, the evolution of the transformed volume is monitored as a 
function of time. Again, we have considered the 100x100 system and introduced a finite 
number (Nq = 20) of nuclei as the initial condition. Therefore, the number density of the 
initial nucleus is n = 20/10000 = 0.002. 

The time evolution of the transformed volume / is plotted as a function of time t in Fig. 3. 
When the effect of the incubation time with to = 10 is included, a better agreement between 
the simulation and theoretical results is attained. This incubation time to = 10 is the time 
necessary for a infinitely small nucleus to become a larger nucleus with the diameter d = 8 in 
our simulation, and is estimated by fitting the theoretical curve (14) to the simulation data. 
Since an infinitely small nucleus cannot grow because it is smaller than the critical nucleus, 
we use the terminology "incubation time" to indicate both the time necessary for a critial 
nucleus to appear and the time necessary for it to grow to be a larger nucleus. 

The KJMA plot of the double logarithm of the volume fraction / is shown as a function 
of log t in Fig. 4, where we ignore the effect of the incubation time and set to = 0. The time 
evolutions for several combinations of the potential parameters r] and e are shown. They do 
not fit the expected straight lines. 
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Fig. 4. Evolution of volume fraction / calculated by cell dynamics simulation for site-saturation 
nucleation as a function of time when t = 0. The double logarithm log (— fn(l — /)) is plotted as 
a function of the logarithm of time logt. All data do not fit the straight lines predicted from the 
KJMA theory (15). 



Figure 5 shows the KJMA plot when the incubation time to is considered. Again, the 
incubation time to is estimated by fitting the theoretical curve (14) to the simulation data. 
Now, the time evolutions for several combinations of the potential parameters r? and e all fit 
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the straight lines with almost the same Avrami exponent n — 2, which is very close to the 
theoretical prediction, as shown in Table II. The results in Figs. 4 and 5 clearly suggest that 
the incubation time to should be carefully taken into account to deduce the Avrami exponent 
n when we analyze experimental as well as simulation data. 
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Fig. 5. Same as Fig. 4 when incubation time to is taken into consideration. All data follows an almost 
straight line and confirm the prediction of the KJMA theory given by eq. (15) with all the Avrami 
exponents n deduced from the straight lines in the figure being close to the theoretically predicted 
n = 2. 



Table II. Avrami exponent n for site-saturation nucleation estimated by cell dynamics simulation for 
various potential parameters n and e. The value theoretically predicted from the KJMA theory is 
n = 2. 



V 


1 


1 


0.4 


e 


0.1 


0.3 


0.1 


n (simulation) 


1.92 


2.04 


2.03 



Figure 6 shows the evolution of the morphology of the two-dimensional system for the 
site-saturation nucleation when r\ = 0.4 and e = 0.1. We observe the almost isotropic growth 
of every nucleus of the stable phase. At the time step ~100, almost all cells are transformed 
into the stable phase. 

3.2.2 Continuous nucleation 

In the continuous nucleation, a new nuclear embryo is continuously introduced. The KJMA 
theory of continuous nucleation gives the analytical expression for the volume fraction / of 
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Fig. 6. Typical evolution pattern of site-saturation nucleation calculated by cell dynamics simulation 
when i] = 0.4 and e = 0.1. The white area indicates the stable phase. 

the growing stable phase. In two dimensions, it leads to 14 

/ = l-exp(-^(t + i ) 3 ), (16) 

where n is the steady nucleation rate per unit area and v is the growth rate of the single 
nucleus discussed in the previous section. 
From eq. (16), we have 

log (- ln(l - /)) = 3 log (t + t ) + constant. (17) 

Therefore, a double logarithmic KJMA plot should give the "Avrami exponent" n = 3 instead 
of n = 2 of the site-saturation nucleation. 

We have also simulated the continuous nucleation using the cell dynamics method. In 
our simulation, a constant nucleation rate h is achieved by introducing a new nucleus every 
1/n time step (nucleation time). At each nucleation time step, a position within the two- 
dimensional area is randomly selected. If the position is already occupied by the stable phase, 
no new nucleus is placed. If the position is not occupied by the stable phase, a new nucleus is 
placed and allowed to grow there. In this simulation, we have used a larger 200x200 system 
to avoid the finite-size effect as much as possible. The steady nucleation rate n = 0.1/40000 
is used. Therefore, a single nucleus is produced at every 10 time steps in the area 200x200. 

The time evolution of the transformed volume / is plotted as a function of time t in Fig. 
7 as the double logarithmic KJMA plot. The time evolutions for several combinations of the 
potential parameters rj and e show almost the same straight line with the Avrami exponent 
very close to the theoretically predicted n = 3, as shown in Table III. 

Figure 8 shows the evolution of the morphology of the two-dimensional systems for the 
continuous nucleation when 77 = 0.4 and e = 0.1. Because a nucleus is continuously produced, 
almost all cells are occupied at the later stage, and then the production of a nucleus stops. The 
situation becomes closer to the site-saturation nucleation. The Avrami exponent n smaller 
than the theoretically predicted n = 3 for the continuous nucleation but closer to the site 
saturation nucleation n = 2 is expected. This finite-size effect is one of the reasons why the 
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Fig. 7. Evolution of volume fraction / calculated by cell dynamics simulation for continuous nucle- 
ation as a function of time. The double logarithmic KJMA plot is used. All curves fit the straight 
lines predicted from the KJMA theory (17). The Avrami exponents deduced from the straight 
lines in the figure are all close to the theoretically predicted n = 3, as summarized in Table III. 

Table III. Avrami exponent n for continuous nucleation estimated by cell dynamics simulation for 
various potential parameters ij and e. The value theoretically predicted from the KJMA theory is 
n = 3. 



V 


1.0 


1.0 


0.4 


e 


0.3 


0.1 


0.1 


n (simulation) 


2.35 


2.60 


2.73 



Avrami exponent estimated from the simulation data in Table III is always smaller than the 
theoretically predicted n = 3. 

There is also a problem of incubation time in continuous nucleation. In our simulation, 
we have introduced a fairly large nucleus, which is sufficiently large to grow continuously. 
Thus, the same problem of time origin or incubation time to as that for the site-saturation 
nucleation could occur. Since a nucleus is continuously produced, we could not incorporate 
the effect of incubation time in a reasonable manner in our analysis. 

In continuous nucleation, we found that the KJMA theory correctly describes the overall 
behavior of the time dependence of the transformed volume fraction /. However, there is 
a small inflection in the slope of the KJMA plot at around log(— ln(l — /)) rs 0.0, which 
can be explained by the impingement in which the growing circular grains collide with each 
other. The volume fraction / when the impingement starts to occur is roughly estimated to 
be / = 7rr 2 /4r 2 « 0.79 using the ratio of the area of a circular grain with the radius r to that 
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Fig. 8. Typical evolution pattern of continuous nucleation simulation when n = 0.4 and e = 0.1. 

of a square with the side length 2r. Thus, the inflection of the KJMA straight line is expected 
at 

log(- ln(l - /)) = log(- ln(0.21)) « 0.19, (18) 

which is very close to the point where the inflection is actually observed in Fig. 7. The same 
effect was discussed by Jou and Lusk. 14 The other factors affecting the Avrami exponent are 
discussed in the next subsection. 



3.3 Discussion 

Experimentally, a reasonably linear KJMA behavior was observed in the recrystallization 
of some metals and in the crystallization of metallic glasses. 9 However, a considerable variation 
in Avrami exponent n defined by 

log (- ln(l - /)) = nlog (t + t ) + constant (19) 

was observed from the electrical resistivity and differential scanning calorimetry (DSC) data. 
Following the argument of Christian, 29 Price 9 proposed a formula for the Avrami exponent 

n = a + b(l-q), (20) 

where a is the nucleation component; a = for the site saturation and a = 1 for the continuous 
nucleation. b defines the dimensionality of the growth (b = 3 for a three-dimensional problem 
and 6=2 for our two-dimensional problem). The exponent q includes contributions from 
various types of power-law reaction. 

For example, in our cell dynamics model based on the Landau- type free energy (9), the 
driving force of phase transformation comes from the undercooling defined by eq. (10) that 
leads to the linear time-dependent growth of a circular nucleus with the constant front velocity 
v given by eq. (12). However, as more materials are transformed into the stable phase, the 
driving force decreases somehow and the front velocity v is expected to decelerate from eq. 
(12). Thus, it is reasonable to assume a power-law decay of the front velocity 

v oc t' q , (21) 
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which exactly gives eq. (20). Hence, the Avrami exponent becomes smaller than the KJMA 
predicted n = a + 6, as observed in our numerical simulations for the continuous nucleation. 

There are many other factors affecting the exponent n. Possible reasons for the nonideal 
exponent and even the nonlinear growth in the KJMA plot include the nonrandomness of 
the nucleation site and the preferential nucleation, for example, at the grain boundary, 30 the 
effect of the time dependence of the nucleation rate and so forth. The net result of these 
effects is a negative deviation from the KJMA linear plot, 9 which leads to again the smaller 
exponent n in accordance to the many experimental results and our simulation. In our cell 
dynamics method, these effects can be easily included by changing the probability of selecting 
the nucleation rate from cell to cell. Hesselbarth and Gobel 19 have included such effects in 
their cellular automaton and could successfully explain the deviation of experimental data 
from KJMA predicted data. 

There is also a problem of two-stage crystallization. 9 ' 31 In some alloys, the KJMA plot 
shows an inflection in which the exponent n changes markedly from a large value at an early 
stage to a small value at a later stage. However, our simulation data shown in Fig. 7 shows 
the opposite trend; the exponent n is large at the later stage. This phenomena is explained 
by assuming that the early stage corresponds to the continuous nucleation and that the later 
stage corresponds to the site saturation because of the exhaustion of the nucleation site. 31 
Recent theoretical model calculation supports this two-stage nucleation model. 32 ' 33 

Our cell dynamics method could easily incorporate such a two-stage transformation by 
assuming that the continuous nucleation terminates at a certain stage. Then, the growth 
process is continuous nucleation with the exponent n = 3 in the early stage, but it becomes 
site saturation with the exponent n = 2 in the later stage. In our cell dynamics method, it is 
not necessary to assume a discrete lattice 32 ' 33 and is easier to incorporate various modifications 
to KJMA kinetics. Using our cell dynamics method, a more quantitative study is feasible in 
the future. 

4. Conclusion 

In this study, we used a cell dynamics method to test the validity of the Kolmogorov- 
Johnson-Mehl- Avrami (KJMA) kinetic theory of phase transformation. First, we used this 
method to study the growth of a single circular nucleus and found that the nucleus grows 
with a constant front velocity in accordance to the analytical solution. 17 Next, we used the 
cell dynamics method to simulate the growth of an ensemble of nuclei under the conditions 
of both the site saturation and continuous nucleation. We found a nearly linear behavior of 
the KJMA plot with the Avrami exponent close to the KJMA predicted value. Finally, we 
suggested several extensions of the cell dynamics method to study various contributions that 
may lead to the nonlinear KJMA behavior or nonideal Avrami exponent. 

The results obtained in this study are summarized as follows: 
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• The cell dynamics method with a realistic free energy can succesfully simulate the steady 
growth of a single nucleus and confirm the prediction of Chan 17 based on the time- 
dependent Ginzburg-Landau (TDGL) equation. 

• It can also simulate the growth of multiple nuclei and confirm the time evolution of the 
volume fraction of the transformed material predicted from the KJMA kinetic theory 1 
and numerical simulation using TDGL. 14 

• Therefore, the cell dynamics method can be used to simulate more complex scenarios of 
nucleation and growth. 

• Our simulation indicates that the incubation time should be carefully taken into account 
when we deduce the Avrami exponent from the experimental and simulation data. 

The cell dynamics method is similar to the time-dependent Ginzbug-Landau or Cahn- 
Hilliard model based on the free energy functional. In contrast to the conventional cellular 
automaton approach to the phase transformation, 6 ' 19,32 ' 33 no phenomenological energy that 
induces phase transformation is necessary. Therefore, the cell dynamics method is numerically 
efficient as a cellular automaton, yet it keeps the direct connection to the equilibrium phase 
diagram. This cell dynamics method can be used to test various scenarios of nucleation and 
growth in a unified manner. 
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